Uncertainty in the Fluctuations of the Price of Stocks 
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We report on a study of the Tehran Price Index (TEPIX) horn 2001 to 2006 as an emerging 
market that has been aflected by several political crises during the recent years, and analyze the 
non-Gaussian probability density function (PDF) of the log returns of the stocks' prices. We show 
that while the average of the index did not fall very much over the time period of the study, its day- 
to-day fluctuations strongly increased due to the crises. Using an approach based on multiplicative 
processes with a detrending procedure, we study the scale-dependence of the non-Gaussian PDFs, 
and show that the temporal dependence of their tails indicates a gradual and systematic increase in 
the probability of the appearance of large increments in the returns on approaching distinct critical 
time scales over which the TEPIX has exhibited maximum uncertainty. 
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I. INTRODUCTION 

In recent years, financial markets have been a focus 
of physicists' attempts for applying the existing knowl- 
edge from statistical mechanics to economic problems [1- 
3] . The markets, though largely varying in the details of 
their trading rules and the traded goods, may be charac- 
terized by some generic features of the time series that 
describe the fluctuations in the prices of various stocks 
and commodities. An important and challenging prob- 
lem is to understand and evaluate risk in the markets, 
which must be done through the analysis of such time 
series. The aim of the analysis is to characterize the sta- 
tistical properties of the time series, with the hope that a 
better understanding of its underlying stochastic dynam- 
ics would provide useful information that can be used for 
creating new models, that are able to reproduce exper- 
imental facts (i.e., the actual recorded prices and their 
fluctuations). 

A considerable amount of data and numerous stud- 
ies indicate the possibility that the financial time se- 
ries may exhibit self-similarity (and/or self- affinity) at 
short time scales which, however, apparently breakdown 
at much longer times. Such features are usually modeled 
in terms of various statistical distributions with trun- 
cated tails. Recent studies indicated, however, that an 
approach based on the Brownian motion [5,6], or other 
more elaborated descriptions, such as those based on the 
Levy and truncated Levy distributions [1], may not be 
suitable for properly describing the statistical features of 
the fluctuations in the stocks' price. Such models have 
been constructed based on the premise that the financial 
time series may be viewed as additive processes that are 



built up over time. There is now increasing evidence that 
an approach based on multiplicative processes might be 
a more fruitful way of pursuing an accurate analysis of 
the financial time series. This approach lends itself in a 
natural way to mult ifr act ality [7] (see below). Such an 
idea was, in fact, suggested some years ago when the in- 
termittency phenomenon in the returns fluctuations was 
observed at different scales, which motivated some efforts 
for establishing a link between analysis of the financial 
time series and other areas of physics, such as turbulence 
[8-11]. We remind the reader that, if pi represents the 
value of a stochastic variable at (time) i, the returns 
are defined by, = ln(pi+i/pi). Nowadays, however, we 
know that there are important differences between the 
two phenomena, such as, for example, the differences be- 
tween their spectra of frequencies. 

Based on the recent efforts for characterization of the 
various stages of the development of markets [12-14], it 
is clear that Tehran stock exchange represents an emerg- 
ing market. It has witnessed considerable activities over 
the past several years, but it is still far from an efficient 
and developed market. Over a two-year period, it lost 
more than 30% of its value (from 13750 units in Septem- 
ber 2004 to 9150 units in August 2006) and, on average, 
the price of the stocks' units has decreased from $0.92 to 
$0.49 (which, percentage-wise, represents an even steeper 
decline than that of the units by which the market has 
lost value), even without considering the rate of infia- 
tion. In addition, over the past six months alone (up to 
the time of writing this paper), the volume and values 
of the traded stocks have decreased by more than 60%. 
Compared with the S&cP 500, Tehran stock exchange is 
still not a completely developed market [13], with its in- 
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FIG. 1: Top: history (2001-2006) of daily deflated closure of 
TEPIX. Bottom: one day log returns of the Tepix index. 



FIG. 2: Continuous formation of increment probability dis- 
tribution function's across scales for, from top to bottom, 
s = 4, 8, 12, 16 and 20 days. The solid lines are the approxi- 
mated PDF based on Casting's equation, the right-hand side 
of Eq (4). 



dex exhibiting stronger non-stationary features. 

In this paper, we provide comprehensive evidence of 
the existence of distinct critical time scales over which 
the Tehran Price Index (TEPIX) has exhibited maximum 
uncertainty. Moreover, at several critical times over the 
past few years, Tehran stock exchange has been affected 
rather strongly by several political crises. These features 
provide a good opportunity to test the method of analy- 
sis suggested by Kiyono et al. [15] for an emerging mar- 
ket. More specifically, by analyzing the temporal evolu- 
tion of the index dynamics, we demonstrate the strongly 
the non-Gaussian behavior of the logarithmic returns of 
the TEPIX and scale-dependent behavior (data collapse) 
of their probability density function (PDF). The critical 
time scales are found to be in the vicinity of large index 
movements, consistent with the high probability of mul- 
tiscale events at the critical points. From the observed 
non-Gaussian behavior of the index, we numerically es- 
timate the unexpectedly high probability of a large price 
change near the critical times. Such estimates are of im- 
portance to risk analysis, as they represent a central issue 
for the understanding of the statistics of price changes. 

The rest of this paper is organized as follows. In the 
next section we present the data that we consider and 
describe how we analyze them. The conclusions are sum- 
marized in Section 3. 
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FIG. 3: The scale-dependence of the fitting parameter of Cas- 
taing's equation vs log s. The inset shows the results for 
a{s) — ^ (('■(i) ~ T^Y')- The results indicate that after s = 4 
days, there is crossover in the behavior of a vs s. 



II. ANALYSIS OF THE DATA 



Figure 1 shows the TEPIX over a period of over 4 i 
years, from December 20, 2001 to August 10, 2006. The 
data had been recorded on each trading day. We show 
in the lower panel of Fig. 1 the one-day log returns, 
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FIG. 4: The local temporal variation of A| ^aya o^sr ^ one-year 
period shows a gradual, systematic increase on approaching 
the critical time scales A-G. 

i.e., rs{t) = \n[p{t + s)/p{t)], where s = 1 day. We 
then analyze the PDF of the detrended log returns over 
different time scales. To remove the trends present in 
{x{t)}, where x{t) = Inp(i), we fit x{t) in each subinter- 
val [1 + s{k — 1), s{k + 1)] of length 2s (where k is the 
index of the subinterval to a linear function of t that rep- 
resents the exponential trend of the original index in the 
corresponding time window. After the detrending pro- 
cedure, we define detrended log returns on a scale s as 
A,p{t) =x*it + s)- x*it), where l + s{k~l) <t< sk, 
and x*{t) is the deviation from the fitting function [7]. 

The scale-invariance properties of a fractal function 
Asp(t) are generally characterized by exponents that 
govern the power-law scaling of the absolute moments of 
its fluctuations, i.e., iTi{q, I) — Kql^i, where, for example, 
one may choose m{q, I) = J^t \^sp{t + I) — Asp{t)\'^ . As 
is well-known, if the exponents are linear in g, a single 
scaling exponent H suffices for characterizing the fractal 
properties with, = qH, in which case Asp{t) is said to 
be monofractal. If, on the other hand, the function is 
not linear in q, the process Asp(t) is said to be multifrac- 
tal. Some well-known monofractal stochastic processes 
are self-similar processes with the following property, 

AxsPit) = A^A,p(i), Vs, A > . (1) 

Widely-used examples of such processes are the fractional 
Brownian motion and the Levy walk. One reason for 
their success is, as it is generally the case in experimental 
time series, that they do not involve any particular scale 
ratio [i.e., there is no constraint on s or As in Eq. (1)]. 

In the same spirit, one can try to build multifractal 
processes that do not involve any particular scale ratio. 
A common approach, originally proposed in the field of 
fully-developed turbulence [8,15-18], has been to describe 
such processes in terms of stochastic equations, in the 
scale domain, describing the cascading process that de- 
termines how the fluctuations evolve when one passes 
from the coarse to flne scales. One can state that the 



fluctuations at scales s and A^ are related (for fixed t) 
through the cascading rule, 

AxsP{t) = WxA,p{t), Vs, A > , (2) 

where In(WA) is a random variable. Let us note that Eq. 
(2) can be viewed as a generalization of Eq. (1) with 
H being stochastic. Since Eq. (2) can be iterated, it 
implicitly forces the random variable W\ to have a log 
infinitely-divisible law [19]. It has been demonstrated by 
Castaing et al. [19] that a non-Gaussian PDF with "fat" 
tails can be modeled by random multiplicative processes. 

Thus, let us assume that the increments in the time 
series are represented by the following multiplicative form 
[7]: 

A,p{t) ^ Cs{t) exp[cu,{t)] , (3) 

where and lOs, assumed to be independent variables, 
are both Gaussian random variables with zero mean and 
variances cr^ and A^, respectively. The PDF of Asp{t) 
has fat tails, depending on the variance of lOs-, and is 
expressed by [19]: 

Ps{Asp) = J Fs (^^^ yGs{\n<Js)d\nas , (4) 

where we have assumed that Fs and Gs are both Gaus- 
sian with zero mean and variance as and As, i.e.. 

Thus, we may investigate the time scale-dependence of 
A^. In this case, the equation for Ps{Asp) is referred 
to as Castaing's equation, whose solution converges to a 
Gaussian when A — > 0. 

The fit of the PDF of TEPIX increments to Castaing's 
equation is indeed almost perfect, especially within ±3 
standard deviations, even for a single record. This is 
demonstrated in Fig. 2. Although Eq. (4) is equivalent 
to that for a log-normal cascade model - originally in- 
troduced to study fully-developed turbulence [19] - it ap- 
proximately describes the non-Gaussian PDFs observed 
not only for turbulence, but also in a wide variety of other 
phenomena, ranging from rate of exchange of foreign cur- 
rencies [8], to heartbeat interval fiuctuations [15,20]. Also 
shown in Fig. 2 is the fit of the data for s ~ 20 days to 
a Gaussian distribution, which clearly fails to represent 
the data. 

For a quantitative comparison, we fit the data (over 
the 4 i years interval) to the above function [Eq. (4)] , as 
illustrated in Figs. 2 and 3, and estimate the variance As 
of Gs- As shown in Fig. 3, the standardized (variance = 
1) PDF of the detrended log returns indicates the exis- 
tence of a scaling law in the behavior of A^ as a function 
of s, rather than logarithmic decay which is character- 
istic of classical cascade processes [17-19,21]. Figure 3 
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indicates that, after s = 4 days, there is a crossover in 
the behavior of as a function of s. For comparison, 
we have also calculated the variance {r{t + T)r{t)) (which 
represents the width of the joint probability distribution) . 
The results are shown in the inset of Fig. 3. Similar to 
As, there is a crossover in the behavior of width. 

In the following, we identify a temporal region of com- 
plete departure from the cascade scenario to an instance 
of the critical- like behavior. We evaluate (in sliding time 
intervals [t — At/2, t + At/2]) the temporal dependence 
of A^ . The local temporal variation of A^^4 ^^^^ over a 
one-year period shows a gradual, systematic increase on 
approaching the critical time scales A-G identified in Fig. 
4. It it beneficial to risk analysis to quantify the non- 
Gaussian nature of (detrended) price fluctuations on a 
relatively short time scale (~ 4 days), and not just the 
volatility at larger time scales [1], which is what is nor- 
mally analyzed. The important point is that large values 
of Aj indicate a high probability of a large price change; 
this probability follows a sharp increase with growing A^. 

The critical points are denoted by A to G in Fig. 4. 
To plot the Fig. 4 we chose a moving window with 
length A< = 150 days. It may be interesting to note 
that these time scales are related to the political devel- 
opments in Iran. There was an increasing trend in the 
price index over the time scale A, caused by privatization 
of Iran's industries. B represents the time period from 
February 21, 2003, when the inspectors of the Interna- 
tional Atomic Energy Agency (IAEA) and its Director- 
General, Dr. Mohammed ElBaradei, travelled to Iran, 
to June 16, 2003, when Dr. ElBaradei reported to the 
IAEA's Board of Governors on what the IAEA had found 
in Iran. G represents the restart by Iran of production 
of centrifuges' parts, used in uranium enrichment (UE), 
on July 31, 2004. D is the time period that included 
the European Union's warning to Iran that it would cut 
off the negotiations on May 11, 2005; Iran's subsequent 
declaration on May 19, 2005 that its UE program is ir- 
reversible, and the election of Iran's new president on 
June 26, 2005. E is the time period over which a new 
director for Tehran stock market was appointed, and the 
economic policies of Iran's new president were declared. 
Finally, F is the time at which the IAEA reported to the 
United Nations Security Council Iran's nuclear dossier 
on February 27, 2006. The time t^, represents the time 
at which Iran's rejection of the IAEA demand for stop- 
ping work on the construction of a heavy-water nuclear 
reactor in Arak was announced on February 13, 2005. It 
can be seen clearly in Fig. 4 that, after that time the 
TEPIX entered a critical period that has continued up 
to now. Moreover, as Fig. 4 indicates, similar to most 
major stock markets around the world, the Tehran stock 
market has responded almost immediately to the polit- 
ical events on the dates indicated. As shown in Fig. 4, 
the trends in the TEPIX are essentially stable up to time 
scale C, but beyond G the average uncertainly increases. 
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FIG. 5: (Color online) Probability distribution function of the 
TEPIX returns before time B and after time A t*. 

To check the changing of the nature of fractal distribu- 
tion of the returns, we plot (in a semi-logarithmic graph) 
the PDFs of the one-day returns before and after the 
critical time t*. The results are presented in Fig. 5. 
Relative to a Gaussian distribution, they exhibit sharp 
peaks, but not long tails. In Table 1, we compare the 
means, standard deviations, skewnesses, and kurtosises 
of the returns time series before and after the time t^, , as 
given in Fig. 1. As Table 1 indicates, the mean value 
of the returns is negative after , but positive before . 
Moreover, the variance after t^ is smaller that its value 
before t*, implying that, on average, the investors have 
lost their investments after but, with smaller risk, had 
gained before t*. 

TABLE I: Comparing the general parameters of the one-day 
return series. 





Mean 


Standard Deviations 


Skewness 


Kurtosis 


t <u 


0.00172 


0.006 


0.84 


16.65 


t > u 


-0.00104 


0.004 


-0.66 


12.46 



III. SUMMARY 

Tehran stock exchange provides a good opportunity to 
test the recently developed analyzing method suggested 
by kiyono et al in an emerging market. We characterized 
the non-Gaussian nature of the detrended log returns of 
the Tehran price index from 2001 to 2006 using a model 
based on multiplicative processes, and found the empiri- 
cal evidence that the temporal dependence of fat tails in 
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the PDF of the detrended log returns shows a gradual, 
systematic increase in the probability of the appearance 
of large increments, on approaching distinct critical time 
scales. The results suggest the importance of the non- 
Gaussian behavior at a time scale (4 days) for risk analy- 
sis. If the same characteristics are observed in other stock 
indices, our approach may be applicable to quantitative 



risk evaluation. 
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